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CH ' 1. Introduction 

The QCD equation of state can be derived from the partition function 

Z{T) = Tre- H / T = Y J 9ne- E "' T . (1) 

In lattice QCD with 2+1 flavours Z(T) has been evaluated by the HotQCD [lj 
and Wuppertal-Budapest [2] collaborations producing different results for 
the trace anomaly at temperatures above T = 200MeV, already beyond 
the hadron-quark-gluon crossover |3j. On the other hand, Quark-Hadron 
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duality at finite temperature requires that for confined states Z should be 
determined from all stable hadron states such as those in the PDG book- 
let [3j. This is the idea behind the Hadron Resonance Gas (HRG), a multi- 
component gas of non-interacting massive stable and point-like particles [5] 
which has historically arbitrated the discrepancies between different lattice 
groups [6HH]- Remarkably, the disagreement still persists beyond the ex- 
pected range of validity of the HRG model (see e.g. Fig. [TJ right). 

The special role played by the HRG does not make it a theorem and 
corrections to it are not completely clear as PDG hadrons are composite, 
have finite size and width. On the lattice, the validity of the HRG has been 
checked in the strong coupling limit and for heavy quarks to lowest orders [9] . 
In the usual large Af c -limit (see Ref. [TO] for a review and references therein) 
where hadrons become stable resonances, T/M = 0(1/N C ), the mesons give 
a finite contribution as their mass and degeneracy are finite whereas baryons 
would provide a vanishing contribution. The half-width rule [TTj applied to 
PDG resonances [3] provides compatible uncertainties with current lattice 
calculations [2]. 

To saturate the partition function, Eq. ([1]) with light or heavy quarks a 
large number of highly excited states is needed so relativistic corrections are 
important. Here, we will use the MIT Bag model [12] and the Relativized 
Quark Model (RQM) of Refs. [131 [H] which treats hadrons as extended 
bound states rather than resonances. 

2. Trace anomaly and light quarks 

The trace anomaly measures departures from scale invariance and reads 

am s lz»- t *.(£), (2) 

after using standard thermodynamics relations for the energy density e = 
E/V and pressure p = —Tlog Z/V. For the HRG model we have 

A(T) - -1 f° dM^± f g* (E k -k.V k E k ) 



where E k = \/k 2 + M 2 and ± corresponds to Fermions/Bosons and 

n(M) = J29a@(M-M a ), (4) 

a 

is the cumulative number ( is the step function). Hagedorn proposed that 
the cumulative number of hadrons in QCD is approximately and asymptot- 
ically given by n(M) = Ae M ' Tn where A is a constant and Tjj is the so 
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. Left: Cumulative number n as a function of the hadron mass M (in MeV) 
u, d and s quarks, computed in the RQM [T3KH] and compared to a fit 



Fig.l 
with 

n(M) = Ae M l Ta . Right: Trace anomaly (e — 3p)/T 4 as a function of temperature 
(in MeV). We compare lattice data for asqtad and p4 [T5] (after temperature down- 
shift of Tq = 15 MeV) and stout [TB] actions, with the HRGM computed with the 
RQM spectrum with u, d and s quarks from Refs. [LUIHF 
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Fig. 2. Left: Cumulative number n as a function of the c-quark mass subtracted 
hadron mass A = M — m c (in MeV) with u, d and s quarks, computed in the 
RQM [Ullli] and compared to a fit n(A) = Ae A/T " . Right: Polyakov loop as a 
function of temperature (in MeV). Lattice data from [T7] for the HISQ/tree action 
and [TB] for the continuum extrapolated stout result. We compare lowest-lying 
charmed hadrons from PDG [1], the RQM spectrum with one b quark and a cut-off 
A < 1700 MeV (red line), and A < 5500 MeV and the MIT bag model (m h -► oo) 
with cut-off A < 5500 MeV is shown as a solid (blue) line [IB] ■ 



called Hagedorn temperature. We show results in Fig. Q] both for n(M) and 
A(T) fitted with A = 0.80 and T H = 260MeV and also showing the good 
performance of the HRG below T = 180MeV when the RQM is used. 
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3. Polyakov loop and heavy quarks 

The Polyakov loop is a purely gluonic operator, which in gluodynamics 
becomes a true order parameter as it signals the breaking of the center 
symmetry and deconfinement. Unlike the trace anomaly, there is lattice 
consensus on this observable [161117] so its analysis may be more credible. 
We have shown that in QCD [18, 19| and in chiral quark models [20] a 
hadronic representation exists and is given by (Aq is the gluon field) 

L T = < tec P e '/o 1/T **o> = IJdA^e-^, (5) 

where the cumulative number reads now 

n(A) = 5> hQ e(A - A a>h ) , (6) 

a 

where g^a are the degeneracies and A^ a = M^a — m^ are the masses of 
hadrons with exactly one heavy quark (the mass of the heavy quark itself 
rrih being subtracted). 

The result with u, d and s quarks, computed in the RQM [13^114] when 
the large but finite charmed quark mass, rrih = m c (using 6-quarks does not 
change much), is taken is presented in Fig. [2j We have checked that results 
are not very sensitive to use bottom quarks instead. A fit n(A) = Ae ' H to 
the total contribution produces A = 0.216, 0.209 and T H = 236, 207MeV for 
single-charmed, bottom hadrons for the range lGeV < A < 1.8GeV. The 
results from PDG and RQM are multiplied by a factor L(T) -)■ e c/T L(T), 
with C = 25 MeV, which comes from an arbitrariness in the renormalization. 
The sum rule has been implemented on the lattice recently |21] , 



4. The non-overlapping condition 

In the quantum virial expansion |22j, the excluded volume corrections 
come from repulsive interactions whereas resonance contributions stem from 
attractive interactions. A good example is tttt scattering where one has 
attractive and resonanting states in the isospin / = 0, 1 corresponding to 
the a and p resonances whereas one has a repulsive core in the 1 = 2 exotic 
channel [23, 23] providing a measure of the finite pion size. In contrast, 
the HRG assumes point-like elementary particles. However, in the narrow 
width limit resonances also have a finite size as they become bound states. 
Clearly, when hadrons overlap, the HRG model becomes invalid since the 
Pauli principle blocks many states allowed by colour neutrality. The non- 
overlapping condition corresponds to the inequality for the Co-volume 

d 3 p gj_ 

(2vr) 3 e E ^P)/T ± i 



oov^wv^v, Em ,J,^-J7t - <i- w 
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Fig. 3. Left: Non-ovelapping condition as a function of temperature. For the hadron 
volume we use Vi = Mi/AB with B = (0.166GeV) 4 for the MIT bag volume (blue) 
and also V% = M 3 /a 3 with y/a = 0.42GeV (red). Right: Cumulative number n(A) 
in the MIT Bag model. We include contributions from Qq, Qqq and Qqqq. 



The hadron size can be estimated from the MIT bag model where one 
has [12] Vi = Mi/(AB). In the RQM [HHHj one might compute the size 
directly from the m.s.r. of the wave functions. A meson model of the form 
M = 2p + ar with p ~ 1/r yields after minimizing V = 47rr 3 /3 ~ M 3 /cr 3 . 
In Fig. [3] we see that for T = 160 — 170MeV hadrons overlap and the HRG 
departs from the lattice QCD results (see Fig. Q] left). 



5. Hagedorn and the bootstrap 

The cumulative numbers computed in the RQM exhibit lower thresholds 
for mesons than baryons but the latter dominate due to the larger multiplic- 
ity of qqq than qq states, and eventually an exponential growth characterized 
by a Hagedorn temperature seems to set in (Figs. [1] and (21) • Due to the finite 
number of degrees of freedom both mesons and baryons have a power-like 
behaviour for large masses M 3> \fo producing a dimensional estimate 
riq q (M) ~ M 6 /cr 3 and n qqq (M) ~ M 12 /a e featuring the available phase 
space. An intriguing issue is under what conditions this exponential growth 
goes on high up in the spectrum as initially speculated by Hagedorn [5]. 
In Fig. [3] (right) we show n(A) = n Qq (A) + riQ qqq (A) + n Qqqq (A) + ... in 
the MIT Bag model including also the exotic tetraquark Qqqq states as in- 
dependent hadronic states. The fit yields Tjj ~ 191MeV, complying with 
the bootstrap mechanism proposed long ago [25j[26]. Since some of the 
tetraquark states are of molecular nature, it is unclear if they should be 
incorporated in the cumulative number. This is related to the completeness 
or redundancy of hadronic states, particularly in the PDG as noted in [11| . 
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6. Conclusions 

The thermodynamical analysis of the hadronic spectrum has an increas- 
ing lack of energy resolution for increasing temperatures and a slowly con- 
verging pattern requiring many excited states. On the other hand, lattice 
calculations become difficult at very low temperatures where the main en- 
ergy gaps are found. While this explains why the HRG model works well 
as function of temperature it is not obvious how to systematically compute 
deviations from this simple limit. 
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